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Background 
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FRC  propulsion 


•  Field  Reversed  Configuration  (FRC) 
is  offshoot  of  fusion  research 

—  Ionization  by  rotating  B-field 
—  Pulsed  inductive  jxB  acceleration 
—  Magnetically  insulated,  plasmoid 
accelerated  downstream 

•  Key  attributes 

—  Very  low  mass  (estimate  ~1-2 
kg/kW) 

—  Efficiency  comparable  to  or  higher 
than  Hall  thrusters  (predicted) 

—  Operates  on  diverse  propellants; 

potential  for  multi-mode  application 
—  Pulsed  operation  provides  near¬ 
constant  efficiency  over  wide  power 
range 
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Rotating  Magnetic  Field  (RMF) 
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Rotating  Magnetic  Field 


RMF  FRC  thruster  design  appears 
well  suited  for  low  power 
operation;  different  set  of 
computational  challenges  than 
theta-pinch  RMF  simulation 


Slough,  J,  Kirtley,  D  and  Weber  T,  Pulsed  Plasmoid  Propulsion:  The  ELF  Thruster ,  IEPC-2009-265,  31th  International  Electric 
Propulsion  Conference,  Ann  Arbor,  Michigan,  September  20-24,  2009. 
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Branch-level  FRC 
development  strategy 


*  In-house 

—  Computational  6.1  program  to  develop  numerical  tools  (this 
program) 

—  Experimental  6.2  program  to  generate  validation  data 

•  Briefly  attempted  theta-pinch  configuration,  failed  to  make 
sufficient  progress  (limited  amount  of  data) 

•  Switched  to  rotating  magnetic  field  (RMF)  configuration 

•  Simultaneous  development  of  PI/PPI  systems 

*  External 

—  Leveraging  SBIRs  to  develop  engineering-level  FRC 
models  for  testing 

—  Testing  thrusters  at  MSNW  LLC  (Redmond,  WA)  and 
AFRL/RQRS  (Edwards  AFB,  CA) 
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Experimental  efforts  (1) 
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Project  Description 
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Justification 


*  Some  computational  models  exist 

—  Neutral  flow  dynamics 

—  Fully-ionized  FRC  translation 

—  Collisional  kinetic  of  charge-exchange 

*  The  right  integration  of  numerical  models  does  not  exist 

—  Partially  ionized  magnetized  plasmas  are  difficult  to 
simulate 

—  Modern  diagnostics  include  internal/external  probes  and 
spectroscopy 

—  Simulation  of  realistic  FRCs  requires  specific  physic 
modules  such  as  coil-plasma  interaction  models 


Fundamental  question  is  how  to  scale  FRCs  to  run  efficiently  at 
low  power  levels  -  difficult/impossible  without  simulation 
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Goals 


•  Pre-ionization  -  need  to  generate  some  amount  of  seed  plasma;  for 
multi-plasmoid  operation,  first  pulse  is  the  most  difficult;  possible  to 
consider  many  different  configurations  (e.g.  liquid  phase  propellants) 

—  Very  specific  to  configuration 

—  Introduction  of  multiphase  and  chemistry  greatly  complicate  physics 

•  Formation  -  go  from  seed  plasma  to  strongly  ionized  closed-field 
plasmoid;  kinetic  effects  expected  in  first  stage;  interaction  between 
coil/plasma  drive  process 

•  Translation  -  plasma  remain  almost  fully  ionized;  interested  primarily  in 
collisional  effects  as  plasma  impact  neutral  background 

—  Fairly  straightforward  to  model  with  existing  fluid  tools  (high  enough 
collisionality  for  Maxwellian  assumption) 


Focus  on  Formation  ->  develop  hierarchy  of  tools  to 
understand  plasmoid  formation,  test  spacecraft  /  integration 

and  design  new  thrusters 
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Objectives 


*  Multi-scale  /  multi-physics  capability  -  emphasis  is  on 
high  order  methods  and  consistent  collision  modules 

*  Formation  and  acceleration  dynamics  -  validation  with 
experimental  test  campaign 

*  FRC  stability  and  turbulence  onset  -  fundamental 
physics  of  current  sheets  in  plasma 

*  Collisional-radiative  characteristics  -  detailed  spectral 
signatures  to  reduce  need  to  perturb  plama 


Well-developed  software  framework  will  support  not  only  FRC 
simulation,  but  also  development  of  future  HET  /  plume  / 
general  partially  ionized,  magnetized  plasma  devices 
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Technical  Challenges 


•  Field/plasma  model  —  Plasmas  can  have  very  complex  phase 
space  configurations  and  self-induced  electromagnetic  forces. 
Fundamental  multiscale  problem  since  electrostatic  forces  are 
relevant  in  this  problem  and  force  very  small  timescales. 

•  Collisional  Physics  —  kinetic  of  ionization  (and  especially 
energy  loss)  are  very  sensitive  to  details  of  electron  energy 
distribution  function 

•  Multi-Scale  Effects  —  diffusive  behavior  in  region  of  separatrix 
is  not  modeled  correctly  in  collisionless  code  -  artificial  diffusion 
makes  computational  phenomena  look  real. 

•  System  Complexity  —  realistic  simulations  need  to  couple 
circuit  models  to  plasma,  represent  three-dimensional  /  end  effects, 
address  radiation  loss  to  impurities,  etc. 
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Approach 


•  Hierarchy  of  Plasma  Models 

—  Continuum:  Hugrass  -  Ideal  MHD  -  Hall/Resistive  MHD  -  MF 
—  Kinetic:  Explicit  -  Implicit  Particle-in-cell  (PIC) 

•  Collisional  Physics  —  Consistent  derivations  from  underlying 
cross-sections  for  both  C-R  and  MF  coupling  terms 

•  Code  acceleration  —  Development  of  physics  modules  in  next- 
generation  simulation  framework  at  AFRL/RQRS;  leveraging  both 
software  (algorithmic)  and  hardware  (GPU)  acceleration 
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Plasma  Description 
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9  Evolve  the  particles  position 
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®  e.g.  Particle-In-Cell  models 
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9  Ensemble  average  of  particles 
distribution,  fv(x,v,t) 

9  Evolve  the  distribution 
function 

9  e.g.  Vlasov-Maxwell  models 
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Kinetic  (non-Maxwellian) 


9  The  Boltzmann  eqn: 
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9  Take  the  0rl\  Vr,  2nd  moments  of  the  Boltzmann  Eqn. 


M, 


j  vu<-^-dv  +  ms  I  v,3+1 


dfs 

Ox 


civ  +  qs  I  vn  (E  +  v  x  B) 


dfs 

Ov 


civ 


m. 


J  o 


dfr 

dr 


c/v 


9  Each  moment  of  the  Boltzmann  eqn  gives  an  equation  for  the  moment 
variable,  and  introduces  the  next  higher  moment  variable 

9  This  process  can  go  on  indefinitely 


Can  implement  Boltzmann  equation  multiple  ways: 

•  PIC  (high  dimensionality  -  3D3V;  high  statistical  noise) 

•  Vlasov  (low  dimensionality  -  2D2V  hero  runs;  smooth  solution) 
Leverage  6.2  program  to  provide  codes  and  expertise 
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Moments  of  the  distribution 


9  Modeling  each  particle  velocity  and  position  is  not  practical. 

9  Instead  an  average  is  performed  to  give  a  statistical  description. 

9  Calculate  the  number  of  particles  per  unit  volume  having  approximately 
the  velocity  v  near  the  position  x  and  at  time  t,  distribution  function 
/(v.  X.  t) 
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No  assumptions  on  shape  of  distribution  function  yet.. ..still  a  general  fluid! 
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Multifluid  (MF)  Plasma  Model 


Fluid  Description  (3  equations  for  each  specie) 
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Further  Simplification 


•  If  there  are  no  neutral  species,  sum  all  charged  species  to  get  a 
single  charged  “fluid”  equation 

•  Connect  with  electromagnetism  through  Faraday’s  law  and  suitable 
version  of  Ohm’s  Law 

(Faraday's  Law)  — -  +  V  x  E  =  0 

Ot 


MHD  variants  of  Ohm's  Law 


Ideal  0  =  E  H-  u  x  B 
Resistive  j  =  cr(E  +  u  X  B) 
Hull  }  =  ^E 

Resistive  Hall  j  —  <r(E  +  u  X  B) 


Magnet: ie  field  frozen  in  flow  (no  diffusion) 
Magnetic  field  can  diffuse  into  plasma 
Incorporates  Hall  currents 
Hall  currents  and  magnetic  field  diffusion 


•  Remove  fluid  equations  altogether  and  simply  represent  plasma  as 
finite  resistivity  affecting  Maxwell’s  equations  -  Hugrass  RMF  model 


Need  all  models  in  same  framework  to  facilitate  simulation  at 
multiply  levels  of  fidelity  can  identify  role  of  missing  physics 
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Detailed  C-R 


*  Require  spectral  signatures  to  interrogate  many  regions  of 
the  plasma  without  significantly  perturbing  plasma 

—  Probes  don’t  function  well  in  RMF  environment 

*  Leveraging  6.2  and  6.1  work  to  access  detailed  Ar  models 

—  Based  on  theoretical  and  experimental  databases  from 
NIST,  LXCAT,  etc. 

—  6.2  work  supports  detailed  model  development  /  6.1  work 
has  supported  intelligent  model  reduction  for  computational 
tractability 

*  Groundwork  has  been  laid  to  develop  self-consistent 
inelastic  transport  coefficients  from  same  cross-section 
databases 
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Progress  Update 
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Numerical  Model  Development 


*  Tackling  both  sides  of  plasma  hierarchy 

—  Low-fidelity  Hugrass  model 

•  Quick  turnaround  so  more  useful  for  interfacing  with 
experimentalists 

•  Potential  low-order  model  for  time-parallel  acceleration 
strategy 

-  High-fidelity  Multifluid  MHD  model 

•  Bringing  SoA  multifluid  capability  to  FRC  modeling 

—  High  order  DG  formulation 

•  Working  on  extension  to  three-fluid 

•  Maintaining  consistency  with  C-R  rates 

*  Computational  machinery  for  high-fidelity  will  be 
used  to  build  MHD  (single  fluid)  system 
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Hugrass  model 


•  Implementation  of  Maxwell’s  equations  in  non-vacuum 


•  Based  on  numerical  model  originally  published  in  1981  by 
Hugrass  (J.  Plasma  Physics,  26,  455-464) 

•  Radial-azimuthal  code  (infinitely  long  cylinder) 


Coupled  equations  for  axial  magnetic  field  and  axial 
component  of  magnetic  vector  potential;  solved  spectrally 
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•  No  self-consistent  plasma  transport  -  can  still  define 
spatial/temporal  dependence  of  resistivity  (77) 
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Vacuum  field  of  RMF  FRC 


Discharge  Tube 

Material:  Quartz 


]  dinner- 6  cm 

wall  thickness  =  2  mm 


Antennas  (A-FI) 

Material:  Copper 
Cross-section:  0.5"  x  0.064" 


Flux  Conserver 

Material:  Copper 
d  inner-  2.756  in 
wall  thickness  0.04  in 


rA  b  =  "KO  cos((oRMpt  +  *P) 
rc,D  =  KQ  sin(coRMF  t) 

rE,F  “  KQ  COS(CORMFt  +  <P) 
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Development  of  MFPM 


Kinetic  LTE,  velocity  moments 

ni,— >0  MHD 
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Ideal  MHD  model  is  valid  when: 

®  High  collisionality,  tu/t  1 

o  Small  Larmor  radius,  ru/L  <C  1 

1  /2 

o  Low  Resistivity,  (^)2  ^  ^  1 

Multi-Fluid  Plasma  Model 


o  Less  computationally  expensive  than  kinetic  models 
®  Multi-fluid  effects  become  relevant  at  small  spacial  and  temporal  scales 
o  Finite  electron  mass  and  speed-of- light  effects  are  included 
e  There  is  charge  separation  is  modeled 
o  Displacement  current  effects  are  resolved  in  the  MFPM 
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Discontinuous  Galerkin  (1) 


„  -i  .  .  £Q  OF  a 

*  Conservation  Laws  are  given  by  —  +  —  =  S 

dt  dx 

®  Multiply  by  basis  functions  and  integrate  over  the  volume 
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That  is  an  ODE  of  the  form 
order  TVD  Runge-Kutta 
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Spatial  order  of  2nd-  16th  is  common. 


=  £r(Q),  solved  using  2"d,  3nl  or  4rh 
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Discontinuous  Galerkin  (2) 


•  Riemann  problems  are  solved  at  each  interface  to  compute  fluxes 

•  The  source  of  dissipation  /  dispersion  depends  on  the  Riemann 
solver 

•  Variables  are  allowed  to  be  discontinuous  at  the  cell  interfaces 

Advantages 

-  Method  is  conservative 

-  Compact  stencil  for  flux  evaluation  (easily  parallelizable) 

-  Provides  good  flux-source  coupling  (no  source  splitting  is  needed) 

-  Resolves  fast  oscillations  (e.g.  plasma  frequency)  and  large  gradients 
(e.g.  shocks) 

Disadvantages 

-  Explicit  Runge-Kutta  Discontinuous  Galerkin  CFL  limit  is  <  1/(2p-1) 
where  high  order  schemes  have  higher  p 
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Multiscale  issues 


Characteristic  speeds: 
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Availability 


Structured  DG: 


9  WARPX  code  developed  at 
Washington 

9  Validated  and  with  proven 
results 

9  Can  be  used  to  verify  the 
unstructured  case 


Unstructured  DG: 


9  Currently  being  developed  at 
AFRL 

9  For  complex 

geometries/experiments 

9  Will  include  different  physics 
modules 
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Verification 


Forward  Facing  Step 
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Backward  Facing  Step 
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Implementation  pieces 


o  Unstructured  grid  readers 

a  Unstructured  data  structures 

o  PETCs  library  for  data 
manipulation  and  time 
integration 

o 

q  XML  input  file  reader 

q  Quadrature  rules  integration 

9  Different  Riemann  solvers  for 
numerical  flux  (e.g.  Roe, 

HLL,  HLLE) 

a  Flux  Limiters  for  shock 
capturing 

a  Initial/boundary  conditions 
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RMF  startup  (1) 


0.7  [is 


Exploring  collisionless  limit  for  efficient  FRC  formation 
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RMF  startup  (2) 


Electron  mass  density 
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Next  Steps 


*  Numerical 

—  Develop/run  continuum  plasma  hierarchy 
—  Run  kinetic  plasma  codes 
—  Integrate  with  collisional  physics 
—  Code  acceleration 

•  Physics 

—  FRC  formation  validation 

•  Magnetic  reconnection 

•  Instability  modeling 
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Conclusions 
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Project  Highlights 


*  Leverages  6.1  and  6.2  efforts 

—  Explicit  Particle-ln-Cell  (PIC)  and  implicit  PIC  development  for  Hall  Effect 
Thrusters  (6.2,  AFRL/RQRS) 

—  C-R  model  reduction  (6.1,  Luginsland/Marshall) 

—  Algorithm  acceleration  (6.1,  Farhoo) 

*  Integration  of  plasma  modules  into  framework  code  provides 
tools  for  investigation  of  numerous  plasma  phenomena 

—  Adapt  Multifluid  description  to  study  electrostatic  anomalous  electron 
transport  due  to  fluid  instabilities 

—  Use  Multifluid  /  MHD  description  as  preconditioner  for  implicit  PIC 

—  Detailed  C-R  necessary  to  evaluate  utility  of  reduced  models  (such  as 
Quasi-Steady  State  (QSS)  approximation) 


Common  simulation  framework  has  potential  to  greatly  increase 
transition  opportunities  by  decreasing  integration  time 
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Summary 


*  Research  plan  directly  supports  R&D  objectives  for 
advanced  plasma  propulsion  (i.e.  FRCs) 

*  Progress  made  in  various  areas: 

—  Moving  towards  consistent  plasma  hierarchy 
—  Integrating  additional  physics  modules  for  real  engineering  utility 
—  Need  more  algorithm  work  (accuracy/stiffness) 

*  Fundamental  work,  can  be  applied  to  other  areas 

*  Challenging,  exciting  project.... 

*  ....but  off  to  a  slow  start 
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